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I. INTRODUCTION 

The connections between cosmology and particle 
physics through dark datter (DM) have recently received 
special attention for defining the physics program at the 
TeV frontier. We foresee that the combination of data 
from satellites, direct DM searches, hadron and lepton 
colliders will provide a major breakthrough in our un- 
derstanding of the nature of dark matter and its inter- 
actions in the early Universe. These expectations are 
supported by the fact that there are several extensions 
of the Standard Model (SM), which include a new, sta- 
ble, weakly-interacting massive particle, which may be 
responsible for the observed relic DM in the Universe. 
This particle should become accessible to particle collid- 
ers operating at the TeV energy frontier, as well as to the 
next generations of direct DM search experiments. The 
LHC collider will be first in providing data to address 
the question of whether one of these scenarios is indeed 
realised in nature. If this is the case, it will also gather 
some quantitative information to be related to the relic 
DM density measured from the cosmic microwave back- 
ground (CMB) spectra However, it is understood 
that the LHC data will not be exhaustive in this respect. 
First, it will not be possible to infer, in a model indepen- 
dent way, the relic density to an accuracy close to that 
already achieved by CMB observations. Furthermore, 
there exist classes of models of new physics which the 
LHC may not be able to disentangle and probe in suffi- 
cient details. It is only with the measurements becoming 
available at an electron positron collider, operating at 
centre-of-mass energies of order of 1 TeV, that we shall 
be able to determine the properties of the DM candidate 
particle and of the other particles participating in its in- 
teractions in the early Universe, with sufficient accuracy 
to predict the DM relic density precisely. With these re- 
sults in hand, the comparison of the data from CMB ex- 
periments, direct DM searches and collider experiments 
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would have striking consequences on our quantitative un- 
derstanding of the nature and distribution of dark matter 
in the Universe. 

In these years preceding LHC operation, Supersym- 
metry has emerged as the best motivated theory of new 
physics beyond the SM. It solves a number of open prob- 
lems intrinsic to the SM and, most important to our dis- 
cussion, the conservation of R-parity introduces the light- 
est neutralino, x?j as a new stable, weakly interacting 
particle. CMB data from the WMAP satellite, and other 
astrophysical data, already set rather stringent bounds 
on the parameters of Supersymmetry, if the lightest neu- 
tralino is responsible for saturating the amount of DM 
observed in the Universe. The recently released, five- 
year WMAP data provide a determination of the dark 
matter density as J] CDM /i 2 = 0.110±0.006 0. 

The potential of the LHC and of an e + e~ linear collider 
operating at 0.5 TeV and 1.0 TeV, such as the Interna- 
tional Linear Collider (ILC), for determining the neu- 
tralino relic density, f2 x , in Supersymmetry has been in- 
vestigated in detail in Q. That study selected a set of 
benchmark points, the so-called LCC points, representa- 
tive of various Supersymmetric scenarios and determined 
the tt x probability density function by a scan of the full 
parameter space of the Minimal Supersymmetric exten- 
sion of the SM (MSSM), by retaining those points com- 
patible with the measurements available at the LHC and 
ILC, within their experimental accuracy. 

In this paper we consider one of the Supersymmet- 
ric scenarios defined in [3], for which the neutralino relic 
density is controlled by its annihilation rate through the 
CP-even heavy Higgs pole \X ~> i which in turn cru- 
cially depends on the value of the mass of the boson, 
Mj^a. We study the accuracy of the measurement of the 
relevant properties of the neutral heavy Higgs boson A : 
its mass, M40, width, T^o and decay branching fractions 
as can be obtained from data collected in high luminos- 
ity e + e _ collisions at centre-of-mass energy of 1 TeV, us- 
ing full simulation of the response of a realistic detector 
model and detailed event reconstruction. 
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FIG. 1: Tranverse energy and thrust distributions for signal 
and background. Generator level distributions are plotted as 
histograms, results of Mokka + Marlin simulation and recon- 
struction are given for the signal process as points with error 
bars. All histograms are normalized to unit area. 



II. e+e~ -> H°A° AT LCC-4 WITH FULL 
SIMULATION 

We adopt the LCC-4 benchmark point of Q , which is 
defined in the reduced paramater space of the constrained 
MSSM by m =380 GeV, mi /2 =420 GeV, tan/3=53, 
A=0, SW^)=+1 and M top =I78 GeV. We use Isasugra 
7 . 69 [4[ to compute the physical particle spectrum 
and we get M A o=419.4 GeV, M x o=169.1 GeV and 
M^ 1 =195.5 GeV. These parameters correspond to a neu- 
tralino relic density of il x h 2 = 0.108, as obtained by using 
themicr oMEGAS 2.0 program The e+e~ -> H°A° -> 
bbbb process at y/s — 1 TeV has already been studied 
for LCC-4 using a parametric simulation Q. We now 
perform a detailed study using Geant-4-based full sim- 
ulation Q of the detector response and reconstruct the 
physics objects using processors developed in the Marlin 
framework Q and extend the analysis to both the bbbb 
and bbr + T~ final states. This study adopts the LDC 
detector concept, which is based on a large continuous 
gaseous tracker, a Time Projection Chamber, surrounded 
by a highly granular SiW calorimeter and complemented 
by a high resolution Si Vertex Tracker. The LDC detec- 
tor concept is discussed in detail elsewhere [9|, the design 
is optimised for achieving excellent parton energy mea- 
surements through the particle flow algorithm, and pre- 
cise extrapolation of particle tracks to their production 
point. Both of these features are important to this analy- 
sis, which aims at suppressing backgrounds by exploiting 
the signature 4-6 and 2-b + 2-r final states of the signal, 
and requires good determinaton of energy and direction 
of hadronic jets to attain an optimal resolution on di-jet 
invariant mass. 

Signal events have been generated with Pythia 
6.205 17] + Isasugra 7.69, including beamstrahlung 
effects 18]. At = 1 TeV, the effective e + e~ -> H°A° 



production cross section, accounting for beamstrahlung 
and initial state radiation, is 1.4 fb, BR(A° — > bb) 
= BR{H° -> bb) = 0.87 and BR(A° -> r+r") = 
BR(H° — > t + t~) = 0.13. The main particle pair pro- 
duction backgrounds, Z°Z°, W + W~ and tt, have been 
generated using Pythia. Their cross sections, computed 
using CompHep 4.4.0 (l9|. are 0.17 pb,_3.0 pb and 0.19 pb 
respectively. The inclusive bbbb and 6&t + t~ production, 
after subtracting the contribution of the Z°Z° channel 
and requiring 200 GeV < M bb < 600 GeV, have cross 
sections of 0.63 fb and 0.28 fb respectively. These pro- 
cesses have been generated at parton level using CompHep 
and then hadronised with Pythia. We assume to op- 
erate the linear collider at \/s=l TeV for a total inte- 
grated luminosity of 2 ab _1 , which corresponds to 5 years 
(1 yr = 10 7 s) of operation for a nominal luminosity of 
4 x 10 34 cm" 2 s _1 . 

A loose event preselection based on event variables has 
been applied after generation. Selected signal and back- 
ground events have been passed through the full LDC 
simulation using the Mokka 06-03 program [l(|, an ILC- 
specific implementation of Geant-4. Data are persisted 
using lcio [111 collections and used as input for the sub- 
sequent reconstruction in Marlin. 

Pattern recognition and track fit are performed first us- 
ing Monte Carlo truth information ( "MC truth patrec" ) 
and, for signal events, also genuine full pattern recogni- 
tion ("full patrec"), using the FullLDCTracking package 
based on DELPHI experiment software [12]. The per- 
formances of these two approaches are compared. The 
Pandora v02-00 package is used for particle flow [l3T ]. 
Jet clustering is performed using the DURHAM algo- 
rithm [13] • The jet energy resolution has been stud- 
ied using a simulated sample of single b jets in the en- 
ergy range from 10 GeV to 210 GeV over a polar angle, 
0.4 < 6 < it/2. We get 5E/E = (0.34 ± 0.02)/ Ve © 
(0.015 ± 0.005), which is consistent with the LDC parti- 
cle flow performance specifications. Jet flavour tagging 
is performed using the LCFIVertex package, which de- 
veloped the original ZVT0P tagger [15| and feeds track 
and vertex topological information into a neural network 
to distinguish between b, c and light quark jets. The 
di-jet mass resolution in the bbbb has been improved by 
performing a constrained kinematic fit. We have ported 
the PUFITC algorithm [l]|, developed for the DELPHI 
experiment at LEP2, into a dedicated Marlin processor. 
The algorithm adjusts the momenta of the jets given by 
Pf = e a pM+bpB+cpc where pp is the fitted momentum, 
Pm is the measured momentum, pb and pc are unit vec- 
tors transverse to pm and to each other, and a, b and c 
are free parameters in the fit. The adjusted momenta sat- 
isfy a set of constraints while minimising the fit \ 2 , given 
by Sj (dj — a ) 2 /cr 2 + 6f/cr 2 + cf/cr 2 , where ao is the ex- 
pected energy loss parameter, a a is the energy spread pa- 
rameter and at,, u c are the transverse momentum spread 
parameters. In this analysis, we impose the following 
constraints: p x = p y — and E ± \p z \ = y/s, where 
the last condition accounts for beamstrahlung along the 
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beam axis, z. 

A. The e+e~ -> H°A° -> 6666 Channel 

First we analyse the fully hadronic final state. This 
provides with characteristics four 6 jet, symmetric events. 
The backgrounds can be significantly suppressed using 
6-tagging, event-shape and kinematic variables. We re- 
quire selected events to fulfill the following criteria: to- 
tal recorded energy in the event E to t > 700 GeV, total 
transverse energy Et >350 GeV, total number of recon- 
structed particles N to t >80, number of charged particles 
iV c ha >30, event thrust <0.9 and Y34 >0.002, where Y34 
is the 3 to 4 jet cross-over value of the jet clustering al- 
gorithm. The distributions of some of these variables are 
shown in Figure [T] for backgrounds and signal, for which 
a comparison of the generator-level and reconstructed 
values is also given. After event selection, particles are 
forced into four jets, which are arranged into two di-jet 
pairs, using the pairing which minimises the difference 
between the di-jet masses, Mjj. The kinematic fit is per- 
formed and a cut applied on the resulting di-jet mass 
difference \Mjji — Mjj%\ <50 GeV to eliminate poorly 
reconstructed events. Both di-jet masses are required to 
satisfy Mjj >200 GeV. The event is required to have 
four 6 jets, where a b jet is determined by the following 
criteria: total jet multiplicity A^ot >1Q 5 charged jet mul- 
tiplicity Ncha >5, and b jet probability, Pb, larger than 
0.5. At the chosen working point, an efficiency for b jets 
of 0.79 is obtained, using "MC truth patrec", with suf- 
ficient rejection of lighter quarks to effectively suppress 
the remaining non-6 backgrounds. By using "full patrec" 
without retraining the neural net, we measure a tagging 
efficiency of 0.72 per jet. 
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FIG. 2: Di-jet invariant mass distribution for signal and back- 
ground events selected by the analysis cuts. Kinematic fit and 
jet flavour tagging have been applied. H°A° events, in which 
the incorrect jet pairing (UP) is chosen, are considered as 
background. 

The di-jet mass for signal H°A° events fulfilling the 



selection cuts has a Gaussian resolution of 23 GeV using 
tracks reconstructed with "MC truth patrec" and 27 GeV 
using tracks from "full patrec" before the kinematic fit. 
After applying the kinematic fit the di-jet mass resolu- 
tions become 13.7 GeV and 13.8 GeV, respectively 

After final selection, the sample of events with di-jet 
masses in the region 200 GeV < Mjj < 550 GeV gives 
a selection efficiency for signal 6666 decays of 0.24±0.01 
using tracks reconstructed with "MC truth patrec" and 
0.17±0.01 using "full patrec". The difference is mostly 
caused by the observed drop in b-tagging efficiency. The 
corresponding acceptance for Z°Z°, W + W~ , it and in- 
clusive bbbb background events is 7 x 10~ 5 , 7 x 10~ 6 , 
8 x 10 -4 and 4 x 10~ 3 , respectively. The resulting mass 
distribution is shown in Figure [21 which has two entries 
per event. The signal is described by the convolution 
of two Breit-Wigner functions with a mass splitting of 
1.4 GeV, as predicted for the LCC-4 parameters, con- 
voluted with a double Gaussian resolution function. The 
background is described by a third-order polynomial with 
coefficients determined on background only events. The 
final fit function consists of a linear combination of the 
signal and background functions with four free parame- 
ters: Ma, Tai and the weights of the signal and back- 
ground functions. We get M A = (419. 7± 1.0) GeV and 
= (14.9±2.9) GeV, where the quoted uncertainties are 
statistical only. This result is remarkably close to that 
obtained in the earlier analysis, based on parametric de- 
tector simulation. Using "full patrec" the uncertainties 
on the A boson mass and width increase to 1.3 GeV and 
3.4 GeV, respectively. 



B. The e+e" -> H°A° -> 66r+r^ Channel 

The mixed decay mode 66t + t~ can be isolated by tag- 
ging a 66 di-jet, consistent with originating from either a 
H° or a A decay and analysing the remaining particles 
in the event. We require the events to fullfi.il the follow- 
ing criteria: E tot >400 GeV, 200 GeV< E T <900 GeV, 
40< N tot <180, 15< N cha <100, event thrust<0.8, event 
sphericity >0.1 and V34 >0.005. The event is forced to 
four jets of which two must be tagged as b jets using 
the same criteria as above but the tighter requirement 
Pb > 0.9. The invariant mass of the bb di-jet must satisfy 
300 GeV< M bb <600 GeV, and that of the two remaining 
jets 250 GeV< Mjj <600 GeV. The angle between the 
two b jets and the angle between the two un-tagged jets 
must satisfy -0.8< cos 8 <0. The number of charged par- 
ticles with energy greater than 5 GeV which are not asso- 
ciated to either of the b jets must not exceed six. Finally, 
r tagging is performed. We have developed an algorithm 
which outputs a linear discriminant variable P T based on 
the jet mass, the impact parameter of the leading track, 
and a variable, Pi sol, which measures the jet energy de- 
posited in an annulus around the jet direction. At least 
one of the two non-6 jets must be tagged as a r jet, where 
a t jet must have less than four energetic charged par- 
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tides and must satisfy P T >0.8. To distinguish between 
signal bbbb and bbr + r~ decays, a discriminant variable 
Pdisc is calculated based on the un-tagged dijet energy, 
the number of energetic charged particles not associated 
to either of the two b jets, and P^ AX , the larger of the 
two tau jet probabilities (see Figure [3|) . The event must 
satisfy Pdisc >0.9. After applying these cuts, the ef- 




FIG. 3: Distribution of the discriminating variable adopted 
for separating bbr + T~ events. 

ficiency for signal bbr + T~ decays is 0.14±0.02, that for 
the tt background is_2_x 10 -4 , for Z°Z° and W + W~ is 
3 x 10~ 8 while for bbbb events is 2 x 10~ 6 . The selection 
criteria yield 87 events of signal with 89 of background, 
corresponding to a relative statistical uncertainty of 0.15 
on the determination of BR(if°, A — > t+t~). 



III. FURTHER CONSTRAINTS ON Sl x 

The constraints on LCC-4 derived from this determina- 
tion of the ^4° mass and width and other supersymmetric 
particle mass measurements at the LHC and a 1 TeV lin- 
ear collider, provide a prediction of the neutralino relic 
density with a relative accuracy of 0.18, within the gen- 
eral MSSM [3]. The main contribution to the remain- 
ing uncertainty comes from the weak constraint which 
data provide to MSSM solutions where Sl x is significantly 
lower than its reference value for LCC-4. A detailed 
study shows that these solutions are all characterised by 
large values of the stau trilinear coupling, A T . In the 
MSSM the f coupling to the H° and A bosons scales as 
^ T fsff an<1 At tan /?+/!, respectively. It has been 
proposed to determine A T through a measurement of the 
branching fraction of A , H° — > fifa [23]. In the fun- 
nel region the main neutralino annihilation mechanism 



is x°X° 



A" -> bb and M A < M f , 



Mr 



A decay into fs allowed by CP symmetry is A — > T\T2 
which is kinematically forbidden for the LCC-4 parame- 
ters. However, at large values of \A T \ 1 the H° — > T\f\ de- 
cay gets a sizeable enhancement of its branching fraction. 
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FIG. 4: H° and A decay branching fractions as a function 
of the stau trilinear coupling A T as predicted by HDECAY. All 
the other MSSM parameters have been kept fixed to those of 
the LCC-4 point. 



In this regime, this channel also contributes to the neu- 
tralino annihilation rate through the x°X° ~ ¥ H° — > ftT% 
process, thus lowering the corresponding relic density, 
as observed in the MSSM scans. At the same time, 
a determination of the branching fraction of the decay 
H — ► fifi, allows us to constrain the stau trilinear cou- 
pling. Figure 0] shows the decay branching fractions of 
the A and H° bosons computed using the HDECAY 2 . 
program [2l[ as a function of the A T parameter. Now, 
due to the same final state, a large H° — ► fifi — > tx°t~x° 
yield can be detected by a standard bbrr analysis, such as 
that discussed in the previous section The present study 
shows that the branching fraction for H°, A — > tt can 
be determined to ± 0.15 and that for A , H° — > bb to ± 
0.07, from which a limit \A T \ < 250 GeV can be derived. 
This constraint suppresses the tail at low values of fl x 
bringing the prediction for the neutralino relic density to 
a relative accuracy of 0.08, which is comparable to the 
current accuracy from the WMAP data. 



IV. CONCLUSIONS 

We have studied the e + e~ — * H°A° process at ^/s = 
1 TeV using on fully simulated and reconstructed events 
for a Supersymmetric benchmark point where the mass 
of the j4° boson is 419 GeV and the relic Dark Matter 
density in the Universe crucially depends on its mass and 
width. We find that the analysis of 2 ab _1 of data should 
probide with relative accuracies of 1.0 GeV and 2.9 GeV 
in the heavy boson masses and widths, respectively. The 
The only branching fractions of the t + t~ decay can be determined 
with a 0.15 relative accuracy. These data, in combination 
with other measurements available at the LHC and a 
e + e~ linear collider, allows to infer the neutralino relic 
density in the Universe with a relative accuracy of 0.08. 
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